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Pseudomonas aeruginosa can cause large and prolonged outbreaks in hospitals. We have sequenced and annotated the genomes 
of two multidrug-resistant P. aeruginosa isolates from the same strain obtained 12 years apart from different patients. Genomic 
analysis provided insight on the genes acquired and lost by P. aeruginosa during its spread. 
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Pseudomonas aeruginosa is a common opportunistic pathogen 
in humans that causes serious infections, especially in neutro- 
penic and critically ill patients. It is intrinsically resistant to many 
antibiotics and can become resistant to all of the antipseudomonal 
agents. Clonal outbreaks of multidrug-resistant P. aeruginosa af- 
fecting hospitals are often described (1,2). The University Hospi- 
tal of Besancon (France) has been hit by a 12-year P. aeruginosa 
outbreak with an international high-risk clone of sequence type 
395 that has affected 250 patients. The antibiotic resistance profile 
of the isolates varied over time (3), along with colony morphol- 
ogy, suggesting an adaptation process that included, but was not 
limited to, the antibiotic pressure. We aimed to decipher the ge- 
netic determinants that drive this evolution during its spread into 
the hospital setting. 

The first isolate of the outbreak, DHS01, was isolated in May 
1997 from the nose of a colonized patient in the surgical intensive 
care unit (SICU). This strain was susceptible to all j3-lactams and 
resistant to fluoroquinolones and aminoglycosides, with the ex- 
ception of amikacin. The later isolate, DHS29, was isolated in 
April 2008 from the urine of another SICU patient. It was resistant 
to all j3-lactams and to fluoroquinolones while susceptible to ami- 
noglycosides. Genomic DNA from the DHS01 and DHS29 strains 
was sequenced on an Illumina GAIIx platform to generate 
22,539,003 and 21,692,025 reads, respectively. The reads were as- 
sembled with Velvet version 1.1 (4) into the draft genomes of 
isolatesDHSOl (6,953,788 bp, 106 contigs, and an N 50 of 159,116) 
and DHS29 (7,108,323 bp, 267 contigs, and anN 50 of 66,667). The 
two genomes had similar G+ C contents of 65.8%. Genome anno- 
tation was performed by the NBCI Prokaryotic Genome Annota- 
tion Pipeline (http://www.ncbi.nlm.nih.gov/genome/annotation 
_prok/) and predicted 6,646 genes, 4 rRNAs, and 61 tRNAs for 
DHS01 and 6,807 genes, 4 rRNAs, and 55 tRNAs for DHS29. 
Comparison of the two genomes with Mummer software (5) re- 
vealed 1,446 single nucleotide polymorphisms (SNPs) and 1,197 
insertions/deletions (indels). We determined the specific regions 



of the two isolates using a cross-mapping of raw reads and assem- 
bled contigs using the Maximum Entropy Method (MEM) algo- 
rithm of the BWA software (6). Seventeen genes were specific to 
DHS01 and 136 to DHS29. A gene encoding a 2 '-aminoglycoside 
nucleotidyltransferase was present in the aminoglycoside- 
resistant DHS01 and absent from the aminoglycoside-susceptible 
DHS29. The DHS29-specific genes were grouped into three clus- 
ters, of which two were similar to the regions pal247-pal292 and 
pal381-pal449 of the reference strain PAOl (7), and included the 
coenzyme B12 biosynthesis gene cluster, the flagellum biosynthe- 
sis cluster/ft, and the quorum-sensing-associated genes lasR, rsaL, 
and lasl. The absence of lasR in DHS01 was consistent with the 
iridescent sheen surface of its colonies on rich agar media (8). We 
show here that during its spread into a hospital setting, an epi- 
demic clone of P. aeruginosa not only accumulates numerous 
SNPs and indels but also acquires or loses gene clusters. 

Nucleotide sequence accession numbers. These whole- 
genome shotgun projects of P. aeruginosa strain DHS01 and 
DHS29 have been deposited at DDBJ/EMBL/GenBank under the 
accession numbers AYNE00000000 and AYON00000000, respec- 
tively. 
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